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Abstract 

In modern surgery, a multitude of minimally intrusive operational techniques are 
used which are based on the punctual heating of target zones of human tissue 
via laser or radio- frequency currents. Traditionally, these processes are modeled by 
the bioheat equation introduced by Pennes, who considers Fourier's theory of heat 
conduction. We present an alternative and more realistic model established by the 
hyperbolic equation of heat transfer. To demonstrate some features and advantages 
of our proposed method, we apply the obtained results to different types of tissue 
heating with high energy fluxes, in particular radiofrequency heating and pulsed 
laser treatment of the cornea to correct refractive errors. Hopefully, the results 
of our approach help to refine surgical interventions in this novel field of medical 
treatment. 
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1 INTRODUCTION 



Presently, modern surgery has at its disposal a great variety of different sur- 
gical techniques to heat biological tissue in a localized and safe way. All 
these techniques are based on specific applicators, i.e. devices which extract 
(cryosurgery) or introduce heat (laser, radiofrequency current, microwave or 
ultrasound treatments). Theoretical heat modeling is a very cheap and fast 
methodology to study the thermal performance of these applicators. In fact, a 
lot of previous work has been conducted to model these procedures by using 
the bioheat equation as governing equation (see [1], and references therein). 
This equation is based on the classical Fourier theory of heat conduction and 
is widely used for modeling the heating of biological tissue. Nevertheless, it has 
been hypothesized that for heat transfers on very small time scales the classical 
model will fail, and an alternative thermal wave theory with a finite thermal 
propagation speed could be more suitable to describe these phenomena. 



For the aforementioned reasons, we are lead to modeling the surgical heating 
of biological tissue by means of the hyperbolic heat transfer equation, and 
then compare these alternative results with those obtained from the standard 
Fourier theory (implying the parabolic heat transfer equation). 



In summary, the main theme of this work is to present in depth the mathemat- 
ical and physical background for a general discussion of heat models related to 
heating of biological tissue by means of energy applicators, and then, in par- 
ticular, point out the physical differences between the classical Fourier theory 
and the hyperbolic wave theory. Moreover, as an illustrative example, we give 
some numerical results of the analytical modeling for radiofrequency heating 
(RFH) and for laser heating, each applied to the cornea in order to correct 
refractive errors. Both surgical techniques have in common that they may 
involve high heat transfers during very short exposure times. 



This paper is organized as follows. In Sec. [2] we will set up the mathematical 
groundwork for heat models related to RFH and laser heating with a special 
emphasis on parabolic and hyperbolic heat models. Then, in Sec. Elwe present 
the results for typical laser and RFH interventions in order to substantiate 
the physical discrepancies between both heat models in concrete examples. 
We will also comment on the domain of applicability of both theories and 
their main features. Finally, Sec. H] will present the conclusions and give an 
outlook on interesting future work. 
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2 HYPERBOLIC VERSUS PARABOLIC HEAT MODELS 



The most fundamental relation to model heat transfer is the generally valid, 
i.e. model-independent, equation for thermal current conservation with a given 
internal heat source 5'(x, t): 

k dT 

V-q(x,t) + - — (x,t) = 5(x,t). (1) 



Here, q(x, t) denotes the thermal flux and T(x, t) is the temperature at point 
X G -D in the domain D C at time t G M+. As usual, thermal conductivity 
is denoted by > and diffusivity by a = /c/pc, where pc is the volumetric 
heat capacity, being p the density and c the specific heat of the material under 
consideration. Fig. [1] gives a schematic view of thermal current conservation. 

Fourier's law of heat conduction [2] constitutes the foundation of all classi- 
cal heat models. It proposes that the heat flux q(x, t) is proportional to the 
negative of the temperature gradient (see Fig. [2]): 

q(x,t) = -/c Vr(x,t). (2) 



However, in this relation temporal changes instantaneously affect heat flux 
and temperature gradient. As an immediate consequence any perturbations 
in classical heat models are propagated with infinite speed. 

Assuming that thermal current conservation ([1]) and Fourier's law ([2]) hold, 
directly yields the classical, i.e. parabolic heat transfer equation {PHTE) with 
heat sources: 

-Ar(x,t) + i^(x,t) = ^^(x,t). (3) 



The standard bioheat equation was introduced by Pennes [3] and is based on 
the parabolic heat transfer model, identifying the following explicit contribu- 
tions for the heat sources in a biological system 

5(x, t) = 5,(x, t) + 5p(x, t) + 5^(x, t), (4) 



where the subindex s denotes a surgical heat source {e.g. laser or radiofre- 
quency treatment), p refers to blood perfusion, and m to any source related 
to metabolic activity. 

With the advancement of modern surgery, medical treatments involve pro- 
gressively smaller time scales and higher energy fluxes. For example, for ra- 
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diofrequency and laser surgery PHTE models could become inappropriate. 
The reason is that on small time scales t G [0,r], with a sufficiently small 
value for r > 0, thermal equilibrium of an extended physical system simply 
can not be reached. This contradicts the classical model, where all elements 
of the thermodynamic system interact with no delay through infinite thermal 
speed. It is to be expected that such theories implying infinite speed of per- 
turbations in the heated medium and more realistic wave theories with finite 
speeds differ considerably in their predictions. 

A modified version of Fourier's law ([2]) uses a non-vanishing relaxation time 
T in the dissipative process. This parameter may be interpreted as the fi- 
nite time necessary for the dissipative flow to relax to its steady thermody- 
namic value. The simplest generalization for finite speeds leads to Cattaneo- 
Vernotte's equation [IfH] : 

q(x,t + r) = -A;VT(x,t). (5) 

Note that the classical Fourier law is still included as a special limiting case, 
taking the limit of zero relaxation time for the heat fluxes. In general, how- 
ever, there is an inertial or retardation term which delays any changes in the 
temperature gradient to be transferred to the heat flux. This can also be seen 
for small relaxation intervals r, which yields the expression 

q(x, t) + r^(x, t) = -k VT(x, t). (6) 

Although it appears at first sight that Eq. ([6]) is an expansion which only holds 
approximately for small relaxation intervals r, by considering entropy argu- 
ments Ozi§ik and Tzou have shown that Eq. is also valid on a macroscopic 
scale mj. 

By combining the heat transfer model ([6]) with current conservation ([1]), one 
directly obtains the hyperbolic heat transfer equation (HHTE) with heat 
sources: 

For zero relaxation the hyperbolic heat equation obviously reduces to the clas- 
sical model. An essential feature of Eq. ([7]) is the second-order time derivative 
with coefficient r > 0, which gives rise to wave properties of the equation, and 
therefore implies finite thermal propagation speeds. This can easily be seen by 
considering forwardly propagating waves which provide solutions for Eq. ([7]) 
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of the form 



T 



■wave 



8 



where k is the usual wave propagation vector and uj the corresponding fre- 
quency. This approach is somehow analogous to the analysis of the telegraph 
equation, which describes the propagation of electromagnetic waves in con- 
ducting media [7I8] . Note that in the present case, however, a term propor- 
tional to T(x, t) is absent on the left-hand side of Eq. ([7]). 

By following the procedure of Ref. and substituting Eq. ([8]) into Eq. ([7]) for 
sourceless heat transfer, then requiring that the real-numbered factor multi- 
plying with A/ has to be zero, one directly obtains for the thermal propagation 
speed 



From this it becomes evident that for vanishing relaxation times r and con- 
stant diffusivity a > 0, the propagation speed tends to infinity. As indicated 
before, but seen from a different perspective, in this limiting case the hyper- 
bolic model reduces again to the standard Fourier theory. 

One further distinctive feature of the HHTE model is observed by representing 
the hyperbolic heat flux equation (jS]) in integral form [H] 



Note that Eq. ffTOj) shows that q(x, t) depends on the full history of the tem- 
perature gradient within the time interval [0,t], something which is entirely 
absent in the much simpler expression for the Fourier flux in Eq. ([2]). 

Exhibiting these features, the hyperbolic heat transfer model with its finite 
thermal propagation appears to provide an ideal framework to describe surgi- 
cal interventions such as laser ablation (a procedure named as laser thermo- 
keratoplasty, or short LTK) and RFH of the cornea {CK, conductive kera- 
toplasty), which use very short and high-energetic pulsations (typical time 
periods are 200 fis and 50 fis, respectively). Therefore, in the following sec- 
tion, we will focus on applications of heat models in surgical techniques, in 
particular RFH and laser surgery applied to the cornea, and thereby develop 
a realistic model for the heating of biological tissue by employing the HHTE 
as previously derived in Eq. ([7]). 




r > 0. 



(9) 




(10) 
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3 SURGICAL APPLICATIONS IN HEAT MODELS 



There exist numerous processes in which great amounts of heat are apphed to 
materials in very short exposure times. This section will focus on the model- 
ing of heating biological tissue as it occurs during medical intervention with 
surgical radiofrequency or laser devices to the cornea. 

Laser heating is a process which implies tissue heating caused by absorbing 
the optical energy of a high-energy laser beam [LTK). On the other hand, 
radiofrequency heating (RFH) is ultimately based on Ohm's law, which states 
that currents flowing through a resistor generate heat. It is essentially ionic 
motion in the tissue which will provoke biological heating aimed at producing 
any medical effects. Both are surgical techniques which may not only involve 
small time scales but also high energy fluxes, which makes them an important 
group of applications in which differences between parabolic and hyperbolic 
models could have great effects. 

As we have previously mentioned in Sec. [21 these type of processes represent 
non-equilibrium processes on small time scales {i.e. compared to the relaxation 
time r), since the system requires considerably more time to reach the equi- 
librium state after the initial thermal energy input. For this reason, PHTE 
models, with their infinite thermal propagation speed, may not provide an 
appropriate description for the underlying physical structure, and it may be 
necessary to rely on HHTE models instead. 

3. 1 Laser Heating 

Among the many surgical procedures in which laser heating is employed, we 
focus our attention on corneal laser heating, also referred to as laser thermo- 
keratoplasty (LTK) 0. 

In the case of laser heating, we consider a theoretical model consisting of 
a semi-infinite fragment of homogeneous isotropic biological tissue in which 
the laser beam falls on the entire tissue surface. Fig. [3] depicts the chosen 
model geometry, in which we consider a one-dimensional model with the x- 
axis parallel to the direction of the incident laser beam. Above the surface, for 
X < 0, free thermal convection is supposed to cool the tissue surface. 

As we want to study the problem from the point of view of the HHTE model, 
we have to start off from the governing equation provided by Eq. ((Tj). For 
this one-dimensional laser heating model, we have to solve the corresponding 
heat transfer equation for penetration depth x including the appropriate heat- 
source term and boundary conditions. 
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Neglecting blood perfusion and metabolic activity in Eq. (jl]), and hence Sp = 
Sjn = 0, the heat source S{x,t) will only have contributions from Ss{x,t). 
This surgical source should be obtained from the Beer-Lambert law, which 
empirically states that for radiation the intensity decreases exponentially with 
penetration depth. Including a factor with the temporal dependence to model 
an energy pulse of duration At, this yields 



S{x, t) = {1- R)bEo e-'"' H{t) - H{t - At) 



:ii) 



where R denotes the dimensionless Fresnel surface reflectance, b is the absorp- 
tion coefficient, and Eq is the incident energy flux at the tissue surface [9]. 
As usual, H{t) denotes the Heaviside function. After combining Eqs. ([7]) and 
fill I) , the final form of the governing equation is 



dx 
[l-R) 



, , 1 (dT , , d'^T , ; 
+ (x,t) + r^(x,t) 



k 



H{t) - H{t -At)+T {6{t) - 6{t - At)) 



(12) 



where 6{t) is the Dirac delta function. 

The corresponding initial boundary conditions are given by 

dT 



T(x, 0) =To, 
lim T(x, t) = Tq 



dt 



(x,0) = 



dT h f dT ^ 

_(0.0 = - r-(0.«)+T(0.0-T 



Vx > 
Vt > 

Vt > 



(13) 



where Tq is the initial temperature and is the ambient temperature. Fur- 
thermore, h > is the thermal convection constant at the interface tissue 
given by Newton's law of cooling, which states that the convective fiux of an 
object is proportional to the difference between its own temperature Tq and 
the ambient temperature Ta- 

qiO,t) = h[Ta-TiO,t)]. (14) 



Since at the interface tissue with x = the heat flux is given by the Newton's 
law of cooling, the last condition of Eq. (fT3!) is obtained by imposing this law 
onto the relation for the hyperbohc heat flux Eq. (fTOj) (see also [6]). 

The explicit analytical solution of this problem has been obtained in Ref. [10] 
and is mainly based on the use of Laplace transforms. It also contains for com- 
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thermal properties 


density 
p [kgm~^] 


conductivity 
k [Wm-^K-i] 


diffusivity 


convection 

h [Wm-^K-i] 


1060 


0.556 


1.3695 • lO"'^ 


20 




optical properties 




Presnel reflectance 
R 


absorption coefficient 

h [m-i] 


0.024 


2000 



Table 1 

Thermal and optical properties of the human cornea. 

parison the fully analytical PHTE solution as an application to the thermo- 
keratoplasty technique. For the numerical estimates, well-established thermal 
and optical properties of the cornea were used. All essential physical constants 
for the model are summarized in Tab. 13.11 

Figure H] shows the numerical estimates for temperatures at various tissue 
penetration depths obtained as a function of time [10]. This figure represents 
the hyperbolic (solid line) and parabolic (dashed line) temperature evolution 
at the four locations x = 0.01,0.1,0.5 and 1 mm for the thermokeratoplasty 
technique. 

We can mainly find two differences between both solutions: 

Firstly, according to the parabolic model, for the values x = 0.01 and 0.1 mm 
the maximum temperature has already been reached and is steadily decreas- 
ing. For X = 0.5 and 1 mm, the maximum temperature has not yet been 
reached and is initially increasing. For the HHTE predictions, we observe 
that for all penetration depths in some part of the plot (or the entire plot) 
the temperature is increasing. The great differences in temperatures between 
both models (overall at positions closer to the surface x = 0.01 and 0.1 mm) 
is even more important when the energy application is pulsed, since in every 
new pulse the difference between the initial temperature could become greater 
and greater. 

Secondly, we can confirm that the behavior of both solutions is very different. 
At values x = 0.01 and 0.1 mm, we notice a delay in the abrupt temperature 
drop associated with switching off the laser beam. This delay is related to 
the finite thermal propagation speed in the HHTE model. Using in Eq. (Q 
the numerical value for diffusivity {viz. Tab. 13.11) and for the relaxation time 
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r = 10 s predicts 

mm , , 

t; = 0.117 (15) 



for the speed of these steps, which exactly agrees with the known finite speed of 
the thermal wave in the cornea pH] . This fact demonstrates the wave character 
of the HHTE model. For depths x = 0.5 and 1 mm, however, these drops 
are not observed, because they occur outside the range of the time interval 
considered. 



3.2 Radiofrequency Heating 



Radiofrequency heating (RFH) is a surgical procedure broadly employed in 
many clinical areas such as the elimination of cardiac arrhythmias, the de- 
struction of tumors, the treatment of gastroesophageal reflux disease and the 
heating of the cornea for refractive surgery [llj. In the remainder of this section 
we will focus on a model of the corneal RFH treatment. 

The schematic diagram of the model geometry is shown in Fig. [51 Here, a 
spherical electrode of radius rg is completely embedded and in close contact 
with the biological tissue, which has infinite dimension. In this case, we can 
also use a one-dimensional model (r being the spatial variable), since the 
geometry under consideration displays radial symmetry. 

Again, assuming Sp = Sm = in the bioheat equation (j4]), the governing 
equation results from the combination of Eq. ([7]) and the explicit form of 
the (surgical) heat source S{r, t) which in the case of RFH is a product of a 
radial [I2] and a temporal part 

Pr 

S{r,t) = S{r)H{t) = ^^Hit), (16) 



so that the total applied power outside of the electrode is 

oo 

P= J S{r)dV = j dVl j r'^S{r)dr, (17) 

shell r>ro ro 



where Q is the usual solid angle subtended by a surface. Thus, the power (flux) 
per unit solid angle, dP/dfl, satisfies the usual 1/r^ intensity law. Note that 
H{t) is included to model a non-pulsed source. As a result, for the corneal 
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RFH the governing equation and the boundary conditions are: 

dT 

T{r,0)=To, _(r-,0) = Vr > tq 

hmT(r,t)=To Vt > (18) 

ot-^ T at PqCqTot or 

The last condition in Eqs. f|T8|) contains the density and specific heat of the 
electrode (where r < ro), given by po = 21500 kgm~^ and cq = 132 Jkg~^K~^, 
respectively. This condition is based on the fact that in the interior and on 
the surface of the spherical electrode it must hold 



dT 

V • q(r,t) +PoCo-^(r,t) 



= 0, (19) 

r<ro 



which just states that the heat fiux q spreading into the biological tissue is 
due to the change of the thermal energy content of the electrode. Applying 
Gauss' law to Eq. f|T9l) and integrating over the entire volume r < tq yields 

/47rr^ dT 
q(r, t) ■ Br rl dQ + pQCo^ — {r, t) = 0, (20) 

where we have assumed that the electrode has relatively small radius and 
sufficiently high thermal conductivity, which makes the electrode act like a 
punctual heat source. In the derivation of Eq. (!20|) . this justifies taking pqCq 
to be constant and pulling dT/dt out of the volume integral. 

Finally, by substituting the expression for the hyperbolic heat fiux Eq. (11 01) 
into Eq. (120]) and using VT • e,. = dT/dr, where is the usual radial unit 
vector, one readily obtains the last condition for the electrode in Eqs. ([11 



The analytical solution of this problem and the corresponding PHTE solu- 
tion have been obtained in Ref. [13] as an application for corneal RFH. We 
have visualized the data generated from Ref. [\A\ in Fig. [6] in order to dis- 
cuss the main characteristics of the model. Figure represents the hyperbolic 
(solid lines) and parabolic (dashed lines) temperature profiles along the radial 
axis for the times t = 10, 30 and 50 ms. Similar to the laser-heating case, for 
small times and locations near the electrode surface the temperatures from 
the HHTE are greater than from the PHTE. This fact also becomes impor- 
tant when instead of a non-pulsed RFH application a pulsed power is applied. 
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Moreover, the presence of pronounced peaks in Fig. [6] reveals the wave char- 
acter of the HHTE model and its prediction of a finite heat conduction speed 



4 CONCLUSIONS AND OUTLOOK 

In this paper, we have outlined the fundamentals and main differences between 
the parabolic and hyperbolic model of heat conduction. These differences en- 
courage the use of the HHTE approach in processes in which great amounts of 
heat are transferred to any material in very short times. Specifically, laser and 
radiofrequency heating are two surgical techniques of this type of processes. 
Through the application of HHTE to these surgical techniques, we have shown 
the main characteristics of this model and its differences with PHTE predic- 
tions. 

At the moment, we are working on the analytical modeling of the hyperbolic 
bioheat equation including a source term for blood perfusion {viz. Eq. (jl])). 
This will be less important for non-perfused tissues (such as the cornea), but 
relevant for well-perfused organs such as e.g. the liver. 

On the other hand, we are studying the implications of the hyperbolic heat 
equation for the case of pulsed RF applications. Energy pulses for surgical pro- 
cedures are being employed in such different areas as conductive keratoplasty 
(CK) and RFH to destroy tumors. 

However, we are well aware of the restrictions of purely analytical models, 
mainly due to the lack of mathematical tractability. For this reason, we are 
planning to develop alternative theoretical models based on numerical tech- 
niques (such as finite elements) in order to use the hyperbolic bioheat equation 
in models with more complicated geometry (more realistic electrode and tissue 
geometries with irregular boundaries) and also varying thermal characteristics. 
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Fig. 1. Schematic view of thermal current conservation. The heat source ^(x, t) in a 
physical system {D € simply-connected domain) is related to both, temperature 
change and divergence of the heat flux q. The coefficient k/a gives the volumetric 
heat capacity in terms of thermal conductivity and diffusivity. 




U thermal energy 

Fig. 2. Fourier's law in differential form, shown for heat transfer through an in- 
finitesimal thin plate with thickness dx and cross-sectional surface area A. The 
corresponding flux, thermal energy per time and area, is proportional to the tem- 
perature gradient dT/dx. The sign indicates that heat flows from T + dT to dT. 
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laser beam 
I I I I 



X = 



free thermal 
convection 



X 



' biological tissue 



Fig. 3. A schematic view of the model geometry for laser ablation. The laser beam 
is applied perpendicular to the tissue and defines the x-axis. The tissue is assumed 
to take the shape of the half-plane x > with free thermal convection above. 




Fig. 4. Temperature profiles for laser heating of the cornea for an initial temperature 
To = 35° C and ambient temperature = 20° C. The laser pulse has a duration of 
At = 200 //s with an incident energy flux Eq = 5 ■ 10^ W/m^. The represented 
curves for the HHTE (solid lines) and Fourier (dashed lines) model correspond to 
penetration depths x = 0.01,0.1,0.5, 1 mm. The relaxation time is r = 10 s. 
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Fig. 5. A graphical sketch of the model geometry for radiofrequency heating, show- 
ing tlie intersection for a spherical electrode. The tip of the active electrode is 
represented by a sphere of radius Vq, and the biological tissue extends to infinity. 
Because of radial symmetry, the problem consists of finding T{r, t), imposing specific 
boundary conditions. 




50 60 70 80 90 100 110 120 130 140 150 160 170 180 



radial axis [(xm] 

Fig. 6. Temperature distributions for the cornea along the radial axis for three 
different measured times t = 10, 30, and 50 ms of a non-pulsed RF application 
with total power P = 30 mW and duration At = 600 ms. The initial temperature 
is To = 35° C and the electrode radius ro = 45 iim. The relaxation time is taken 
r = 0.1s. 
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